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ABSTRACT 

Using data drawn from the Sloan Digital Sky Survey (SDSS) and the SDSS~II Supernova Survey, 
we study the local environments of confirmed type la supernovae (SNe la) in the nearby Universe. 
At 0.05 < z < 0.15, we find that SN la events in blue, star-forming galaxies occur preferentially in 
regions of lower galaxy density relative to galaxies of like stellar mass and star-formation rate, while 
SNe la in nearby red galaxies show no significant environment dependence within the measurement 
uncertainties. Even though our samples of SNe in red hosts are relatively small in number, tests on 
simulated galaxy samples suggest that the observed distribution of environments for red SN la hosts is 
in poor agreement with a cluster type la rate strongly elevated relative to the field rate. Finally, after 
considering the impact of galaxy morphology, stellar age, stellar metallicity, and other relevant galaxy 
properties, we conclude that the observed correlation between the SN la rate and environment in the 
star-forming galaxy population is likely driven by a gas-phase metallicity effect, such that prompt 
type la supernovae occur more often or are more luminous in metal-poor systems. 
Subject headings: supernovaeigeneral, galaxies:statistics, galaxies: abundances, galaxies:stellar content, 
large-scale structure of universe 



1. INTRODUCTION 

Type la supernovae (SNe la) are thought to be dis- 
tinct from other types of supernovae (that is, result- 
ing from a different progenitor population), as they are 
found in galaxies spanning a broad range of properties. 
While, type II, lb, and Ic SNe are only found in star- 
forming galaxies, indicating that they are the evolution- 
ary product of massive stars, SNe la are found in both 
star-forming and quiescent systems, suggesting that they 
are someh ow connected to the evo l ution of less-massive 
stars Ce.g.-IOemler fc Tinslevlll979l: Ivan den Berghl [l990l: 



Idella Valle k Livid 119941 : iCaPDellaro et al.1 I1999D . In 

particular, it is widely accepted that the progenitors 
of type la supernovae are carbon-oxygen white dwarfs 
(WDs), which have accreted mass up to the Chan- 
drasekhar limit ()Chandrasekhaij | 193lD. perhaps via de- 
position from a binary companion ( Whelan fc Ibenlfl973t 
iHan fc Podsiadlowskil [2004) . Lending support to this 
picture, th eoretical model s of exploding W Ds (e.g., 
ikasen fc Plewa 2005; Kasen fc WooslevI l2007t ) are able 
to reproduce the properties of SN la spectra, including 
the lack of hydrogen features and the presence of strong 
silicon absorption, which define the type la classification. 

Recent observations of the type la SN rate in local and 
intermediate-redshift galaxies have cast some doubt on 
— or at minimum confused — the theoretical paradigm 
just discussed. If SNe la are simply the product of 
old stellar populations (i.e., an evolutionary outcome 
of low-mass stars), then the SN la rate should depend 
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strongly on stellar mass, while being independent of the 
level of on-going star-formation activity (at fixed stel- 
lar mass) . While the type la rate is found to depend on 
stellar mass, it is also found to be greater per unit stellar 
mass in galaxies with higher specific star-for mation rates 
(jMannucci et al.ll2005HSullivan et al.ll2006b[ ). These cur- 
rent observations of the super nova la rate support a 
revised theoretical mo del (e . g., IScannapieco fc Bildstenl 
120051 iMannucci eFall l200d iPritchet et all l2008l ). p~ 
haps one employing a two-component progenitor distri- 
bution with a "prompt" component correlated with star- 
formation activity and a "delayed" component correlated 
with stellar mass (i.e., the underlying older stellar pop- 
ulation). Still, despite uncertainty in the nature of the 
two observed components of the type la rate as well as 
difficulties in directly observing the progenitors of type 
la events, current observations remain consistent with 
a progenitor population comprised entirely of carbon- 
oxygen WDs, allowing for an increasingly bro ad distri- 
bution of delay times (e.g.. iGreggio et al.ll2008l ). 

Studying the environments of supernova — i.e., their 
place in the hierarchy of large-scale structure — may 
prove useful in shedding light on the nature of the two 
components of the type la rate or on the nature of the 
underlying progenitor population or potentially even re- 
vealing dependencies of the type la rate or luminosity 
on galaxy properties that are correlated with environ- 
ment. Type la supernovae are often used as cosmological 
probes, allowing distances to be measured out to interme- 
diate redshifts and providing constrain ts on cosmological 
parameters such as Hq, VIm, and JIa (jRiess et al.lll998t 
[Perlmuttcr et al. 1999). Understanding any relationship 
between host galaxy properties and type la luminosity 
is critical for minimizing systematic effects that might 
bias studies of distant type la supernovae and limit our 
abilities to constrain cosmological models. Along these 
lines, environment can be utilized as a proxy for other 
galaxy properties that are otherwise difficult to constrain 
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observationally (e.g., stellar metallicity or age). For in- 
stance, a galaxy of the same stellar mass as another, but 
located in a higher-density environment, will generally 
have formed earlier in the history of the Universe and 
will contain older stellar po pulations as a result (e.g., 
ISchave et al1l2003l : iDave et al. 20061 

Understanding if the type la rate varies with the lo- 
cal galaxy environment may provide new insights into 
the true nature of the progenitor population, while also 
shedding light on many aspects of galaxy formation and 
evolution. For instance, using environment as a proxy for 
formation time could be helpful in interpretin g the role of 
an evolving initial stellar mass function (e.g.. lDavell2008l : 
Ivan Dokku m 2008) in establishing the potential progen- 
itor population. 

More directly, constraining the correlation between the 
SN rate (of both type la and II events) and galaxy envi- 
ronment is also critical for understanding various details 
of the chemical enrichment and star-formation histories 
of galaxies. For exa mple, super novac gover n the pro- 
ducti on of metals (Wooslev fc W eaver 1995; iSato et al.l 
|2007|) . with types la and II dominating iron and oxygen 
production, respectively. Due to their large gravitational 
potential wells, galaxy clusters are potentially the only 
systems to have retained all of the metals produced by 
stars. For this reason, studies of the metal content (e.g., 
Fe abundance) in the intracluster medium (ICM) pro- 
vide interesti ng windows on the cosmic s t ar-formation 
histo ry (e.g., iMatteucci &: Vettolanil [l988t ICalura et al.l 
l2007t ). Given the role of supernovae in dispersing met- 
als, a critical part of this picture is to understand how 
the S N rate may vary with e nvironment in addition to 
time (jSivanandam et al.l[200i l. 

Furthermore, energy injected into the interstellar 
medium (ISM) by supernovae could play an impor- 
tant role in influencing t he formation of galactic disks 
(jScannapieco et ahl I2008D . Variation in the super- 
nova rate with environment could thus be a factor 
(though likely very weak) in the establishment of the 
morphol ogy-density relation (e.g.. iDavis" fc Gelled [19761 : 
iDressleii |l980). Finally, feedback from supernova- 
generated winds is thought to be an e ssential part 
of shaping the mass-metallicity relation (|Dekel fc Silkl 
ll986l:ICoig[l99l . causing a downturn in the relative en- 
richment at low stellar masses by ejecting metals into the 
intcrgalactic medium (IGM) . Understanding the correla- 
tion between supernova rates and environment will help 
reveal whether supernova feedback is also responsible for 
driving the scatter in this fundamental relation, causing 
the observed correlation between metallicity and envi- 
ronment at fixed ste llar mass and star-formation rate 
(jCooper et al.ll2008a[) . 

In this paper, we utilize public data froni the Sloan 
Digital Sky Survey (SDSS, IVork et al.\\2W^ and fr om 
the SDSS-II Supernova Survey (jFrieman et al.l 120081 ) to 
examine the environments of local type la SNe. In §2, 
we discuss the data samples employed along with our 
measurements of galaxy environments. Our main results 
regarding the environments of type la SNe are presented 
in §3, with comparison to related work, analysis of po- 
tential selection effects, and further discussion in §4. Fi- 
nally, in §5, we summarize our conclusions. Through- 
out this paper, we assume a flat ACDM cosmology with 
f^m = 0.3, Q\ = 0.7, w = — 1, and a Hubble parameter 



of Ho — 100 km s ^ Mpc ^, unless otherwise noted. 

2. THE DATA SETS 
2.1. The Supernova Sample 

O ver the past three yea r s, the SDSS-II Su pernova Sur- 
vey (jFrieman et all l2008t ISako et all 12008 ') has repeat- 
edly scanned a 300 square degree region around the ce- 
lestial equator in the southern Galactic hemisphere in 
search of supernovae. The resulting data set, which 
includes more than 800 SNe, is exceptional in its size 
and uniformity, significantly increasing the number of 
known supernovae at z < 0.3 with well-sampled and 
well-calibrated light curves. Only the Supernova Legacy 
Survey fSNLS. lAstier et al.l [20061 ). which primarily tar- 
gets supernovae at higher redshifts {z > 0.3), has discov- 
ered a comparable number of type la events. 

In this work, we utilize the > 500 confirmed SNe la 
from the SDSS-II Supernova Survey.^ We limit the sam- 
ple to a redshift range of 0.05 < z < 0.15 and to SNe for 
which a spectrum was obtained, which yields a total sam- 
ple size of 163 type la SNe (SN-A sample). We exclude 
type II supernovae from our analysis due to their limited 
numbers; there is only one type II supernova for every 

10 type la event in the SDSS-II Supernovae Survey cat- 
alog. W hile type II e vents are more numerous than SNe 
la (e.g., iMaod [20081 ) . they are intrinsically fainter and 
thus less-likely to be detected in a magnitude-limited 
survey. 

A key aspect of the SDSS-II Supernova Survey is 
the general spatial uniformity of the number of imag- 
ing epochs. That is, there is little spatial dependence 
to the on-sky cade nce of the survey (se e Figure 1 of 
iFrieman et al[ 120081 and lSako et"aI[[2008D . This trans- 
lates to a survey that is equally sensitive to supernovae in 
all environments, which therefore allows for direct com- 
parison of the environments of SN hosts to the environ- 
ments of the global galaxy population. 

2.2. The Galaxy Sample 

With spectra and multi-band photometry for more 
than 500, OOP galaxies, the SPS S Data Release 6 (DR6, 
lAdelman-McCarthv et all l2008l ) enables the local den- 
sity of galaxies (which we refer to throughout this pa- 
per as "environment") at z < 0.2 to be measured 
over approximately one quarter of the sky, including 
nearly the entire area surveyed by the SDSS-II Super- 
nova Survey. We select a sample of 553, 188 galaxies 
from the SDSS DR6, as contained in the New York 
University Val ue-Added Galaxy Catalog (NYU-VAGC, 
iBlanton et al.|[2005 b). This sample is limited to the red- 
shift range 0.01 < z < 0.3 and to SDSS fiber plates for 
which the redshift success rate for targets in the main 
spectroscopic survey is 80 per cent or greater. 

For each galaxy, the SDSS/NYU-VAGC database pro- 
vides precise information about position on the plane of 
the sky (<0.1 arcsecond rms pe r coordinate) and alon g 
the line of sight {Az ~ 0.0001. [Abazaiian et ai]|2004D . 
enabling the local environment to be accurately charac- 
terized (see §2.4|) . In addition, rest-frame colors and ab- 
solute magnitudes are computed usin g the KCORRECT 
iiT-correction code (version v4_l_4) of lBlanton fc RoweisI 

^ Downloaded from |http:/ /sdssdp47.fnal.gov/sdsssn/sdsssn.html[ 
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(|2007l see also iBlanton et alj l2003"al) . The rest-frame 
quantities for the SDSS sample are derived from the ap- 
parent ugriz model magnitudes in the SDSS DR6, where 
all SDSS magni tudes within this pa per are calibrated to 
the AB system (Oke & Gunh'1983'). 



As shown by man y previous studies at low and interme- 



diate redshift (e.g., 


Strateva et al.ll2001l; iBell et al.ll2004l: 


iWillmer et alJ 12006 


) and as illustrated in Figure [Tl the 



distribution of galaxies in color-magnitude space is bi- 
modal, with a tight red sequence and a more diffuse blue 
cloud. To divide the SDSS galaxy sample into these two 
broad classes, we use the following magnitude-dependent 
cut: 



.g-r^ -0.02667- Af^ + 0.11333. (1) 

This division in rest-frame g — r color is shown in Fig. [1] 
as the dashed red line. 

Stellar masses are computed for each galaxy in the 
SDSS sample, again using the KCORRECT package with 
templat e spectral energy distribut ions (SEDs) based on 
those of iBruzual &: CharlotI (|2003D . To estimate stellar 
masses, the best-fitting SED given the observed ugriz 
photometry and spectroscopic redshift is used to directly 
comp ute the stellar m ass-to-light ratio (M,/L), assum- 
mg a iChabrieil (|2003f l initial mass function. 

Due to the large wavelength coverage of the SDSS spec- 
tra (3800-9200A), star-formation rates (SFRs) are able 
to be estimated for nearly every galaxy in the spectro- 
scopic catalog based on the measured Ha (A6563A) line 
emission. For each galaxy in the spectroscopic sample, 
emission-line fluxes and equivalent widths (EWs) are 
measured by fitting for and subtracting the stellar con- 
tinuum, as detailed bv lYan et al.l ()2006[ ). To derive the 
line luminosity and correct for aperture effects related to 
the finite size of the SDSS fibers, we estimate the total 
Ha luminosity by combining measurements of the Ha 
EW with the isT-corrected broad-band absolute magni- 
tudes (i.e., assuming that the ratio of the line emission to 
the broad-band fiux is uniform across the entire galaxy) . 

Star-formation rates are inferred from the measured 
Ha luminosities according to the relation given by 
iKennicuttI (|1998f) : 

V'(Ha) = 7.9 X 10-42 yr"!, (2) 

ergs s ^ 

where L(Ha) is corrected by a factor of 2.8 to account 
for underlying dust attenuation and stellar absorption. 
In a small percentage of cases (< 15%), the Ha flux 
is unable to be measured accurately due to bad pix- 
els within the emission-line or continuum windows. In 
these instances, we infer the SFR from the measured [O 
II] A3727A line l uminosity, correc t ed us ing the empiri- 
cal calibration of iMoustakas et all ()2006f ). Note that in 
computing our star-formation rates, we adopt a Hubble 
parameter, h — 0.7, to match that calibration. 

Our measured SF Rs agree well with those measured 
for SDSS galaxies bv lBrinchmann et al.l ()2004D , who used 
flts of emission-line and co ntinuum properties to stellar 
population models (see also lCharlot et al.ll2002f l: through 
direct comparison with our inferred star-formation rates, 
we find an offs et of 0.3 dex and a scatt er of ~ 0.15 dex 
relative to the iBrinchmann et aT] (120041) measurements, 
with the offset of the IBrinchmann et "all (12004 ) SFRs to 
higher values largely due to differences in dust correc- 



tions. Direct co mparison of our e s timat ed Ha luminosi- 
ties to those of IMoustakas et all (|2006D show excellent 
agreement, with only a small offset corresponding to a 
^ 0.01 dex offset and ^ 0.02 dex scatter in the inferred 
star-formation rates. 

In addition to the SDSS spectroscopic data set, we 
utilize the larger SDSS DR6 photometric catalog, which 
contains uniform, precision photometry for millions of 
sources down to a 5-(t limiting magni tude of r = 22.2 in 
asinh mag nitudes (jLupton et al.l ll999t). For all sources m 
the imaging catalog, we estimate rest-frame colors, ab- 
solute magnitudes, stellar masses, and photometric red- 
shifts using the KCORRECT package. Due to the lack 
of spectral information, we are unable to estimate star- 
formation rates for sources in the imaging catalog. 

2.3. Identifying Host Galaxies 

For each supernova in the SDSS-II Supernova Survey 
sample, we attempt to identify a host galaxy within the 
SDSS DR6 spectroscopic and imaging data sets. When 
identifying hosts in the SDSS spectroscopic galaxy cat- 
alog, a projected, radial window of 25 kpc (phys- 
ical) is employed to distinguish potential hosts on the 
plane of the sky in conjunction with a velocity window 
of Av — 3000 km s~^ along the line of sight. This 
moderately-large velocity window is adopted to account 
for the relatively-low precision of the redshifts derive d 
from SN la features {Az < 0.005. iFrieman et al.ll200i l. 
If multiple SDSS galaxies fall within the radial window, 
then the galaxy closest in projected distance is taken as 
the host. A host is considered to be unambiguously iden- 
tified, if one (and only one) galaxy is found within the 
search window. 

In the redshift range 0.05 < z < 0.15, a total of 48 
type la supernovae (SN-B sample) are matched to host 
galaxies (45 of the 48 are matched unambiguously) in 
the SDSS spectroscopic catalog. The location of the 
host galaxies in color-magnitude space is given by the 
green points in Figure [T^; the hosts are roughly equally 
divided between the red sequence (26/48) and the blue 
cloud (22/48). We exclude one supernova from the SN-B 
sample due to its occurrence in a bright, nearby QSO (see 
the hexagon point in Fig. [T|) , for which measurements of 
luminosity and star-formation rate are highly uncertain. 

Although the SDSS spectroscopic data set supplies rel- 
atively precise information about the line-of-sight posi- 
tion of many galaxies in the area surveyed by the SDSS- 
II Supernova Survey, the SDSS imaging catalog provides 
a significantly more complete census of the galaxy popu- 
lation due to its much greater depth; the SDSS imaging 
catalog is complete down to r = 22.2, while the SDSS 
spectroscopic sample is magnitude-limited at r < 17.77. 
For this reason, we also search for host galaxies within 
the SDSS DR6 imaging catalog. 

Using a magnitude-limited (r < 19.5) imaging cata- 
log,® we identify potential hosts on the plane of the sky 
within a projected, circular window of 25 kpc (phys- 
ical) in radius. To differentiate between potential hosts 
along the line-of-sight, we compute the redshift for each 
galaxy in the imaging catalog, within the radial window, 

^ We employ a relatively bright magnitude limit, to ensure high- 
precision photometry and photometric redshifts, thereby minimiz- 
ing contamination by other objects along the line— of— sight. 
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Fig. 1. — {Left) We plot the rest— frame g — r versus M,- color— magnitude distribution for SDSS galaxies in the spectroscopic sample 
within the redshift range 0.05 < z < 0.15. Due to the large number of galaxies in the sample, we plot contours (rather than individual 
points) corresponding to 25, 50, 200, 500, 1000, and 2000 galaxies per bin of A(g — r) = 0.05 and AMr = 0.1. The dashed red horizontal 
line shows the division between the red sequence and the blue cloud as given in Equation [T] The green points correspond to the 48 SN la 
host galaxies included in sample SN— B, where the two galaxies within lh~^ comoving Mpc of a survey edge are outlined by red diamonds. 
The cyan hexagon identifies one particular host, which is a low-redshift QSO. The blue circles outlining three of the green points denote the 
three SNe with ambiguous host identifications. The magenta stars (connected by dotted lines to the corresponding blue circle) illustrate 
the location in color-magnitude space of the other potential hosts. In one case, there are three possible hosts in the spectroscopic sample. 
Finally, note that the SN hosts are divided roughly evenly between the red sequence and the blue cloud. (Right) Same as plotted on the 
left, except that the green points show the rest-frame g — r versus Mr color-magnitude distribution for the 60 SN la host galaxies in sample 
SN— C. Again, the two galaxies within 1 comoving Mpc of a survey edge are outlined by red diamonds and the one SN hosted by a 
QSO is marked by the cyan hexagon. 



using the photometric-redshift c ode SDSS.KPHOTO Z 
in KCORRECT (version v4_l_4 iBlmTton et al.l [2003ah . 
Within a velocity window of Ati = 6000 km s~^, we se- 
lect the closest galaxy in projected distance as the host. 
A total of 60 type la supernovae in the redshift range 
0.05 < z < 0.15 (SN-C sample) are matched to host 
galaxies in the SDSS imaging catalog. The distribution 
of these hosts in color-magnitude space is given by the 
green points in Figure [TJd; they are weighted more to- 
wards the red-sequence population relative to the hosts 
of the SN-B sample (36 out of 57 are red). The SN- 
C sample is a superset of the SN-B sample, with all of 
the supernovae in SN-B being matched to the same host 
galaxy. Again, we exclude the one supernova from the 
sample that is found within a local QSO. 

2.4. Measuring the Local Environment 

We consider the "environment" of a galaxy to be de- 
fined by the local mass overdensity, as traced by the lo- 
cal overdensity of galaxies; over quasi-linear regimes, the 
mass density and galaxy density should simply differ by 
a factor of the galaxy bias (Kaiseii.1987). To estimate 
the overdensity of galaxies in the SDSS, we utilize mea- 
surements of the projected fifth-nearest-neighbor sur- 
face density (S5) about each galaxy, where the surface 
density depends on the projected distance to the fifth- 
nearest neighbor, Dp,^, as S5 = 5/(7rZ?p 5). In computing 

S5, a velocity window of ±1500 km s~^ is employed to 
exclude foreground and background galaxies along the 
line of sight. The projected distance to the N*''~nearest 
neighbor provides an accurate estimate of local galaxy 
density o ver a broad and cont inuous range of scales. As 
shown bv lCooper et al.l ()2005l ). it is reasonably robust to 
redshift-space distortions, while also effectively tracing 
the local density in underdense regions. 

To correct for the redshift dependence of the sampling 
rate of the SDSS, each surface density value is divided by 
the median S5 of galaxies at that redshift within a win- 



dow of A2 = 0.02; this converts the S5 values into mea- 
sures of overdensity relative to the median density (given 
by the notation 1 + 5^ herein) and effecti vely accounts 
for re dshift variations in the selection rate (jCooper et al.l 
l2005f ). We restrict our analyses to the redshift range 
0.04 < z < 0.16, avoiding the low- and high-redshift 
tails of the SDSS dN/dz distribution where the varia- 
tions in the survey selection rate are greatest. Finally, 
to minimize the effects of edges and holes in the SDSS 
survey geometry, we exclude all galaxies within 1 
Mpc (comoving) of a survey boundary, reducing our sam- 
ple size to 392, 938 galaxies within the redshift range 
0.04 < z < 0.16. 

For each supernova (those with and without an identi- 
fied host in the SDSS DR6 galaxy catalog), we measure 
the local environment in a manner identical to that fol- 
lowed for the galaxy sample. That is, we measure the 
local surface density of galaxies about the position of 
the supernova, using the positional information (a, S, z) 
from the SDSS-II Supernova Survey and using the SDSS 
DR6 galaxy sample to trace the local environment. After 
excluding SNe near the survey boundary (within 1 
Mpc), we arrive at a final sample (SN-A) of 134 SNe la. 
Note that for the subset of SNe with an identified host 
in the SDSS DR6 spectroscopic catalog (SN-B), higher- 
precision information about the local environment is also 
available by proxy via the host galaxy. A summary of all 
of the supernova samples is provided in Table [TJ 

2.5. Selecting Comparison Samples 

As discussed in [JT] a variety of recent observations have 
shown that the SN la rate (weighted by mass or by lumi- 
nosity) de pends on the properties of the host galaxy. For 
instance, 'Mannucc i et al.l (20051). using the supernova 
catalog of Cappella ro et al. ( )1999l ). showed that SNe la 



are more common in morphologically late-type galaxies 
(Irr and Sbc/d) relative to more bulge-dominated sys- 
tems (E/SO). Similarly, recent work from the Super- 
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nova Legacy Survey (| Sullivan et aLl l2006a[ ) found that 
the SN la rate is greater among galaxies with greater stel- 
lar m ass and among galaxies with higher star-formation 
rates (jSullivan et al . 2006b). 

For several decades, the observed properties of 
galaxies (including star-formation rates, morphol- 
ogy, and rest-frame color) have been known to de- 
pend u pon the local envi r onmen t [e.g.. iDavis fc Gelled 
197a (Postman fc Gelleil fl98l IBaloeh et all Il99a 



Cooper et al.ll2006[ ). In particular, galaxies with more 



massive stellar populat ions tend to favor regions of 
highe r galaxy density (|Hogg et al.l 12004 IZehavi et al.1 
I2OO5I ) , while systems with high star-formation rates typ- 
ically reside in low'-dens i ty environs in the lo cal Uni- 
verse (jGomez et al.ll2003t ICooper et al.|[2008bl) . In ad- 
dition, the relationships between galaxy properties and 
environment depend on redshift, with the color-density 
and morpholog y-density relations growing weaker at 
higher z (e.g iDsg ssler^eLalJ 119971: iSmith et al.l 120051 : 
iCooper et al 1 120071 ) and the SFR-d ensity relation in- 
verting between z y 1 and z ~ (jElbaz et al.l l2007l : 
ICooper et aLll2008bD . 

Given these known correlations between galaxy prop- 
erties and [1] the SN la rate as well as [2] the local galaxy 
density, we extract multiple subsamples from the SDSS 
galaxy sample, selected to match the characteristics of 
the SN la samples, thereby enabling analysis of the local 
environments of SN hosts independent of selection biases 
connected to galaxy type. Since many of the SNe in the 
SN-A sample lack an identified host galaxy in the SDSS 
galaxy catalog, we are unable to select a comparison sam- 
ple that matches properties such as luminosity, stellar 
mass, star-formation rate, etc. Here, we randomly select 
15,000 galaxies from the set of 392,938 SDSS galaxies 
with accurate environment measures, so as to match their 
redshift distribution to that of the SN-A sample. We uti- 
lize a matching radius of Az — 0.01, randomly drawing 
galaxies from the redshift range 0.04 < z < 0.16. In Fig. 

we show the redshift distributions for this subsample 
(SDSS-A) relative to that of the SNe in the SN-A sam- 
ple. 

For the smaller SN-B sample, however, we are able 
to randomly select a comparison sample that matches 
galaxy properties such as luminosity, color, and stel- 
lar mass. From the set of 392,938 SDSS galaxies with 
accurate environment measures, we draw two compari- 
son samples: one matched to the luminosity, color, and 
redshift distributions of the SN hosts in sample SN-B 
(sample SDSS-B) and a second matched to the stellar 
mass, star-formation rate, and redshift distributions of 
the same host galaxies (sample SDSS-F). 

Members of the comparison samples are drawn ran- 
domly from within 3-dimensional radial windows of 
A{g - r)2 + AAf2 + Az"^ < 0.002 and Alog(V')^ + 
Alog(M,)^ -I- Az^ < 0.002, centered on the properties 
of each host. The SDSS-B and SDSS-F samples are 
constructed from 7, 500 independent, random matches. 
Thus, some SDSS galaxies are duplicated in the compar- 
ison samples; however, the large size of the SDSS spec- 
troscopic galaxy catalog ensures that duplication is min- 
imized such that > 80% of the comparison samples are 
unique, with no individual galaxy included more than 8 
times in a given comparison sample. 

Figure [3] shows the relative distributions of rest-frame 
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Fig. 2. — The redshift distributions for the 134 SNe la in sample 
SN-A and for the 15, 000 galaxies in sample SDSS-A. As summa- 
rized in Table[T] the SDSS-A galaxy sample is randomly chosen to 
match the redshift distribution of the supernovae. The red dotted 
and blue dashed-dotted lines show the redshift distributions for 
red and blue subsets of the SDSS-A sample, respectively, following 
the division in rest-frame g — r color as given by Equation [T] Note 
that all histograms are normalized to have equal area. The total 
number of objects per bin for the SDSS-A and SN-A samples is 
shown by the left and right axis scales, respectively. 

color, luminosity, and redshift for the galaxies in the 
SN-B and SDSS-B samples; by design, the distribu- 
tions of these properties are well matched. Given the 
relatively tight relationship between the combination of 
rest-frame optical color and lum inosity with stella r mass 
and star-formation rate (e.g., iKauffmann et all 120031 : 
ICooper et aT1l2008bf l. we find that the distributions of 
stellar masses and star-formation rates for the SN-B and 
SDSS-B samples are also closely matched (see Figure |4]). 
The SDSS-F sample, which is directly matched to the 
stellar mass and star-formation rates of the host galax- 
ies in the SN-B sample, shows similar distributions of 
these galaxy properties. 

Finally, we define a comparison sample (SDSS-C), se- 
lected to match the distribution of rest-frame colors, lu- 
minosities, and redshifts of the SN-C supernova sample. 
Recall that the SN-C sample is selected by matching the 
SDSS-II supernovae to the SDSS DR6 imaging catalog. 
Due to the lack of spectroscopic information for all of the 
hosts, we are unable to estimate accurate star-formation 
rates for the SN-C sample. The details of each of the 
supernova and galaxy samples is given in Table [TJ 

3. RESULTS 

In Figure [5l we show the environment distribution for 
the 134 supernovae in the SN-A sample in comparison to 
that of the SDSS-A galaxy sample. The SNe appear to 
be clearly biased towards lower-density environs relative 
to a sample of random galaxies with the same redshift 
distribution. Moreover, the overdensity distribution for 
the supernovae looks to be skewed towards low overdensi- 
ties even relative to that of the blue galaxies in SDSS-A, 
where the blue subsample is selected following the rest- 
frame color division given by Equation [TJ 

A variety of statistical tests have been developed 
to determine whether two sets of data come from 
the same underlying distribution. We have applied 
two of the most powerful non-parametric tests (i.e., 
those which are independent of Gaussian assump- 
tions) to our dat a: the one-sided Wilc o xon-Mann- 
Whitney U test (|Mann fc WhitnevI 119471: iPress et all 
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TABLE 1 

Descriptions of the SNe Ia and Comparison Galaxy Samples 



Sample A^obj A''edge —cut ^ range Brief Description 



QIV A 

oiN — rt. 




1 f\'i 


1 Q/1 


U.UO <. z \J 


1 

Id 


all OiN 6 ia witii spectra irom tne ouoo— ii oi\ ourvey 


SDSS- 


-A 


15,000 


15,000 


0.04 < ^ < 


16 


galaxies chosen randomly to match the redshift distri- 
bution of the SN-A sample 


SN-B 




48 


45 


0.05 < ^ < 


15 


all SNe la with spectra from the SDSS-II SN Survey 
and with identified host galaxies within ij < 25 
kpc in the SDSS spectroscopic sample, selecting the 
closest host in projected distance in cases of confusion 


SDSS- 


-B 


7,500 


7,500 


0.04 < ^ < 


16 


galaxies chosen randomly to match the luminosity, 
color, and redshift distributions of the SN— B sample 


SN-C 




60 


57 


0.05 <z <0 


15 


all SNe la with spectra from the SDSS-II SN Survey 
and with identified host galaxies within R < 25 
kpc in the SDSS imaging sample 


SDSS- 


-C 


7,500 


7,500 


0.04 < ^ < 


16 


galaxies chosen randomly to match the luminosity, 
color, and redshift distributions of the SN— C sample 


SN-D 




48 


45 


0.05 < ^ < 


15 


all SNe la with spectra from the SDSS-II SN Survey 
and with identified host galaxies within i? < 25 
kpc in the SDSS spectroscopic sample, selecting the 
bluest host in cases of confusion 


SDSS- 


-D 


7,500 


7,500 


0.04 < ^ < 


16 


galaxies chosen randomly to match the luminosity, 
color, and redshift distributions of the SN-D sample 


SN-E 




48 


45 


0.05 < ^ < 


15 


all SNe la with spectra from the SDSS-II SN Survey 
and with identified host galaxies within R < 40 
kpc in the SDSS spectroscopic sample, selecting the 
closest host in projected distance in cases of confusion 


SDSS- 


-E 


7,500 


7,500 


0.04 < ^ < 


16 


galaxies chosen randomly to match the luminosity, 
color, and redshift distributions of the SN— E sample 


SDSS- 


-F 


7,500 


7,500 


0.04 < ^ < 


16 


galaxies chosen randomly to match the stellar mass, 
star-formation rate, and redshift distributions of the 
SN-B sample 



[l99l IWallfc Jenki^ [20031) a nd the Kolmo eorov- 
Smirnov test ( Press et al. nTool I Wall fc Jenkhlil2003l) . 



Note. — We list each supernova and galaxy sample employed in the analysis, detailing the selection cut 
used to define the sample as well as the redshift range covered and the number of objects included before 
(A^obj) S'lid after (Af(,dgo— cut) removing those within 1 comoving Mpc of a survey edge. 

ther sense) less appropriate. Since the test is one-sided, 
possible P-values range from to 0.5. The Wilcoxon- 
Mann- Whitney (WMW) U test is particularly useful for 
small datasets (as we have for our SN samples) due to 
its insensitivity to outlying data points, its avoidance of 
binning, and its high efficiency. 

In addition to the U test, we also employ the two-sided 
Kolmogorov-Smirnov (KS) test. This test quantifies dif- 
ferences between samples using the maximum absolute 
difference between the cumulative distributions of two 
datasets, the K-S statistic. This difference will be small 
for data drawn from the same distribution and large for 
dissimilar data; it is particularly sensitive to differences 
in the "core" of a distribution, but relatively insensitive 
to differences in tails, lending the test high robustness to 
non-Gaussianity. The P-value for the test is then the 
probability of obtaining the observed value of the K-S 
statistic, or a higher one, for the null hypothesis that we 
have two random datasets drawn from the same under- 
lying distribution. If this probability is small, then the 
likelihood is high that the distributions of the data are 
not identical. 

Performing a one-sided WMW U test on the overden- 
sity measures for the SN-A and SDSS-A samples, we 
find that the SN-A environment distribution is skewed 



The result of each test is a P-value: the probability that 
a value of the U or K-S statistic equal to the observed 
value or more extreme would be obtained, if some "null" 
hypothesis holds. Throughout the remainder of this pa- 
per, results with a P-value below 0.05 (corresponding 
closely to 2a for a Gaussian) will be considered to be 
"significant", while P-values below 0.01 will be classi- 
fied as "highly significant" . 

The Wilcoxon-Mann- Whitney U statistic is computed 
by ranking all elements of the two datasets together, and 
then comparing the mean (or total) of the ranks from 
each dataset. Because it relies on ranks, rather than 
observed values, it is highly robust to non-Gaussianity, 
but still has efficiency (as measured by the sample size 
required to reach a given error level) almost as high as 
the classical t test for Gaussian distributions. In partic- 
ular, we apply a one-sided U test, which determines the 
P-value for the null hypothesis that a specified sample 
is not skewed to higher values than the other sample; 
the obvious differences between the environments of the 
blue supernova hosts and blue galaxy samples (e.g., see 
Figure [7]) make a two-sided test (for the null hypoth- 
esis that the samples have the same distribution in ei- 
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Fig. 3. — The rest— frame color (g — r), absolute magnitude 
(Mr), and redshift (z) distributions for the 45 SNe la host galaxies 
in sample SN-B and for the 7, 500 galaxies in sample SDSS-B. All 
histograms are normalized to have equal area. The total number of 
objects per bin for the SDSS-B and SN-B samples is shown by the 
left and right axis scales, respectively. As summarized in Table [T] 
the SDSS-B galaxy sample is randomly chosen to match the color, 
luminosity, and redshift distributions of the supernovae hosts. 

to smaller values of overdensity than that of the SDSS-A 
galaxy sample, with a P-value < 0.01, i.e., there is less 
than a 1% chance that we would observe a difference this 
strong if both samples were drawn from the same parent 
distribution. If we compare to the blue SDSS-A galaxy 
subsample, we still obtain P = 0.02, which is highly sig- 
nificant. If instead we simply compute means and errors 
on the mean for the various overdensity distributions in 
Fig. O we arrive at this same general result (see Table 
[2]). The mean overdensity, < logjQ(l -l-Js) >, for the SNe 
in sample SN-A is notably lower than the mean environ- 
ment of the SDSS-A sample, with the difference signifi- 
cant at a > 3 CT level. Relative to the mean overdensity of 
just the blue galaxies in the SDSS-A sample, the typical 
environment of the supernovae is still underdense, but 
with the difference only significant at a ^ 1 cr level. A 
KS test confirms these results, indicating that the SN-A 
and SDSS-A environment distributions are very unlikely 



to have been drawn from the same parent distribution. 

The differences in the environment distributions for the 
SN-A and SDSS-A samples, however, are potentially 
only driven by differences in the properties of the two 
galaxy samples in concert with the underlying correla- 
tions between galaxy properties (e.g., color, morphology, 
etc.) and environment. For instance, a sample of galax- 
ies chosen randomly from the SDSS DR6 galaxy catalog 
will undoubtedly turn out to be dominated by red galax- 
ies; the red fraction is 60% for the full SDSS sample at 
0.05 < z < 0.15.^ If we assume that the hosts of the SNe 
in the SN-A sample have colors that relatively closely fol- 
low the distribution of the SN-B hosts (roughly evenly 
divided between red and blue), then our SDSS-A sam- 
ple is likely somewhat more skewed towards red galaxies 
than are the hosts of the SN-A sample. Since galaxy 
color is strongly correlated with environ ment, with red 
galaxi es favoring overdense regions (e.g., iBlanton et al.l 
l2005al) . the SDSS-A sample is potentially biased towards 
a sample of galaxies in higher-density regions, compared 
to the SN hosts, just because of this effect. For this rea- 
son, meaningful conclusions cannot be drawn from the 
analysis of the SN-A and SDSS-A samples. 

The SN-B sample, which includes only those SNe with 
identified host galaxies in the SDSS spectroscopic sam- 
ple, allows us to compare the environments of type la 
supernovae to the environments of a sample of galaxies 
with similar properties to those of the SN hosts, thereby 
avoiding the biases that weaken the utility of the larger 
SN-A sample. As discussed in §2.5) the SN-B sample 
has the additional benefit of more precise environment 
estimates for the supernovae (relative to the estimates 
for the SN-A sample), since the host galaxy redshifts 
have roughly 50 times higher precision than the super- 
nova redshifts. 

In Figure m we show the environment distribution for 
the 45 SNe la in the SN-B sample alongside that for 
the SDSS-B sample, which was selected to match the 
color, luminosity, and redshift distributions of the 45 SN 
hosts. As we found for the SN-A sample, the SNe la in 
the SN-B sample appear to be more commonly found in 
lower-density regions relative to both the blue and the 
red galaxies in the SDSS-B sample. Performing a WMW 
U test on the overdensity measures for the SN-B and 
SDSS-B samples confirms that the SN-B environment 
distribution is distinct from that of the SDSS-B galaxy 
sample with a P-value less than 0.04, a significant result. 

By dividing the SN-B sample according to the rest- 
frame color of the host galaxy (again, following Equation 
[1]) , we are also able to directly compare the environment 
distributions for the blue and red galaxies separately. In 
Figure [71 we plot the distributions of overdensity mea- 
sures for these subsamples selected by rest-frame color. 
While the distributions for the red galaxies in SN-B and 
SDSS-B appear to be quite similar, the overdensities for 
the blue SN la hosts are significantly skewed to lower val- 
ues relative to the blue galaxies in the SDSS-B sample. 
Wilcoxon-Mann- Whitney U tests confirm these impres- 
sions, concluding that the overdensity measures for the 
blue hosts are stochastically smaller than those of the 

^ The red fraction increases with redshift in the SDSS sample 
and is actually > 60% at z 0.12, the median redshift of the 
SDSS-A sample (see Fig.EJ. 
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Fig. 4. — The distribution of stellar masses (left) and star-formation rates (right) for the 47 SNe la host galaxies in sample SN-B and 
for the 7, 500 galaxies in sample SDSS-B. Each histogram is normalized to have equal area, with the total number of objects per bin for 
the SDSS-B and SN-B samples indicated by the left and right axis scales, respectively. As summarized in Table [T] the SDSS-B galaxy 
sample is randomly chosen to match the color, luminosity, and redshift distributions of the supernovae hosts; however, the stellar mass and 
star— formation rate distributions are well matched too. This reflects the strong correlations between specific SFR and stellar mass-to-light 
ratio with rest-frame color. 

sity measures for SN-B to those of the SDSS-F sample, 
which is a galaxy sample selected to match the stellar 
masses, SFRs, and redshifts of the SN-B host galaxies, 
rather than their colors, luminosities, and redshifts. As 
shown in Tableland Table [3l we find that the blue SN 
hosts in SN-B are biased to lower-density regions rel- 
ative to the blue galaxies in the SDSS-F sample, with 
a P-value (or significance) similar to that found when 
comparing SN-B to SDSS-B. 

While the environment measures for the SN-C sample 
are less precise than those of SN-B, due to the lack of a 
spectroscopic redshift for each host galaxy, we still detect 
a significant difference between the environment distribu- 
tion of blue hosts in the SN-C sample in comparison to 
like galaxies in the SDSS-C sample. The mean overden- 
sities for the blue subsamples differ at a 3.5cr level, while 
the typical environments of the red subsamples are indis- 
tinguishable within the uncertainties. The WMW U and 
KS tests support the results derived from analyzing the 
mean overdensities, with the environment distribution 
for the blue hosts distinguishable from that of a color-, 
luminosity-, and redshift-matched sample at P < 0.01 
and P < 0.03, respectively. 

Altogether, the primary result of analyzing the large- 
scale environments of samples SN-B and SN-C is that 
blue SN la host galaxies are found to be biased towards 
low-density environments relative to galaxies of like stel- 
lar mass, star-formation rate, and redshift. As discussed 
in more detail in 21 this result can be interpreted as 
evidence for a bias in the rate or luminosity of type la 
events in low-density regions, such that prompt la events 
are more numerous or more luminous in underdense en- 
virons. 

One of the most striking features of the environment 
distribution for the blue hosts in SN-B, as shown in Fig. 
[6l is the complete lack of blue host galaxies in high- 
density regions. All of the supernovae in star-forming 
systems are found in overdensities of log]^Q(H-55) < 0.14, 
while the blue comparison galaxies in the SDSS-B sam- 
ple span the full range of overdensities. To test the sig- 
nificance of this sharp cut-off in overdensity, we draw 
100, 000 random subsets of 20 galaxies each from the blue 
SDSS-B sample. Of these 100, 000 subsamples, wc find 
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Fig. 5. — We plot the distribution of overdensities for the SN- 
A (green dashed line) and SDSS-A (black solid line) samples. In 
addition, we divide the SDSS-A sample according to rest-frame 
g — r color, following Equation [T] and plot the environment dis- 
tributions for the resulting red (red dotted line) and blue (blue 
dashed-dotted line) subsamples. All of the histograms are normal- 
ized to have equal area, with the total number of objects per bin 
for the SDSS-A and SN-A samples indicated by the left and right 
axis scales, respectively. We find that the environment distribution 
for the SNe la in sample SN-A is skewed to lower overdensities rel- 
ative to a random sample of galaxies from the SDSS with a similar 
redshift distribution. 

blue SDSS-B galaxies at a > 99% level, while the distri- 
butions for the red SN-B and SDSS-B galaxy samples 
are indistinguishable from each other. Examining the 
mean overdensities, < log]^Q(l -I- S5) >, for the respec- 
tive samples confirms the result of the WMW U tests. 
As shown in Table O the mean environments for the 
blue samples differ at a > 4 tr level, with the SN hosts 
being typically found in lower-density regions relative to 
galaxies of like color, luminosity, and redshift. Two-sided 
Kolmogorov-Smirnov (KS) tests similarly find that the 
blue hosts in the SN-B sample have an environment dis- 
tribution distinct from that of the blue SDSS-B galaxies 
with P < 0.01, while no significant distinction is found 
when comparing the environments of the red hosts in 
SN-B and the red galaxies in SDSS-B. 

The difference between the environments of the host 
galaxies in the SN-B sample and the environments of like 
galaxies is further supported by comparing the overden- 
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that less than 0.1% display sharp cut-offs in their en- 
vironment distribution such that all 20 galaxies reside 
in overdensities of log;^o(l + 65) < 0.15. This analysis 
further supports the conclusion that the SN la rate or 
luminosity is elevated in low-density environments (rel- 
ative to more overdense environs) among star-forming 
galaxies. 

We test the strength (or depth) of this cut-off in over- 
density by measuring how much of an increase in the 
type la rate in low-density environments is needed to 
make the observed cut-off at log]^Q(l -I- 55) = 0.15 sta- 
tistically likely to occur. To do this, we select the 3244 
random blue galaxies in the SDSS-B sample, which by 
design match the rest-frame color, absolute magnitude, 
and redshift distributions of the blue SN la hosts in 
sample SN-B. From this parent population, we then 
draw 100, 000 independent samples of 20 random galax- 
ies with the likelihood of drawing an object in a low- 
density (log;^Q(l + 65) < 0.15) environment forced to 
be 2x greater than that of a galaxy in a high-density 
(log]^Q(l-|-(55) > 0.15) environment. The division between 
the underdense and overdense regimes is selected to be 
logiQ(l-|-^5) = 0.15, so as to maximize the probability of 
seeing an apparent cut-off. We repeat this exercise while 
varying the degree to which the SN rate increases in the 
underdense regime (e.g., 3x and 4x greater, etc.). 

We find that to have a 1%, 5%, or 32% probability of 
observing a cut-off at logiQ(l -I- ^5) = 0.15 would require 
an increase in the SN la rate in low-density environments 
at the level of >1.5x, 2.5x, and 7x, respectively. Thus, 
a SN la rate ^ 2x higher within star-forming galaxies 
in low-density environments (relative to those in more 
overdense regions) is likely required to produce the ob- 
served cut-off in the environment distribution among the 
blue SN host galaxy population. 

From this statistical analysis, we have determined the 
probability of observing so strong a strong cut-off in 
environment, given a SN rate in high-density environs 
which is R times that in low-density regions, as a func- 
tion of R. From this, we can determine the Bayesian 
equivalent of a 95% confidence interval, a 95% credible 
interval, as this will be the region contai ning 95% of the 
posterior probability Bayes' theorem (.Press et al.lll992l : 
iWall fc JenkinsI I2003D shows that this posterior proba- 
bility will be proportional to the probability of obtaining 
a cut-off as strong as observed (or stronger) for a given 
R — the likelihood — multiplied by the probability dis- 
tribution we would assign to R in the absence of any 
measurements — a prior. Based on this analysis, we 
conclude that there is 95% probability that R < 0.38; 
i.e., that the SN la rate is less than 0.38 times as large 
in the high-density regime (logjQ(l + ^5) > 0.15) as at 
lower densities, assuming a flat prior probability distri- 
bution for R, a conventional zero-information prior for 
a parameter with some characteristic scale (e.g., R = 1). 
If we instead adopt a prior with uniform probability for 
all intervals of log(_R) (i.e., P{R) — 1/i?), we would find 
that there would be only a 5% probability that R > 0.12. 
Finally, even with a prior as extreme as PiR) = R, which 
favors i? > 1, we find that there is 95% probability that 
R < 0.66; an equal SN rate in high- and low-density 
regions (i? = 1) is strongly ruled out by our analysis.* 

* Note that this analysis assumes that the supernovae in high- 
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Fig. 6. — We plot the distribution of overdensities for the SN— B 
(green dashed line) and SDSS— B (black solid line) samples. In ad- 
dition, we divide the SDSS— B sample according to rest-frame g — r 
color, following Equation [T] and plot the environment distributions 
for those red (red dotted line) and blue (blue dashed-dotted line) 
subsamples. All of the histograms are normalized to have equal 
area, with the total number of objects per bin for the SDSS-B 
and SN-B samples indicated by the left and right axis scales, re- 
spectively. We find that the environment distribution for the type 
la SN hosts in the SN-B sample is skewed to lower overdensities 
relative to a randomly— selected sample of galaxies with matching 
luminosity, color, and redshift distributions. 

4. DISCUSSION 

In SjHl we show that the SNe la in blue host galaxies oc- 
cur preferentially in low-density environments relative to 
galaxies of like color and luminosity (or like stellar mass 
and star-formation rate). For the type la events in red 
hosts, however, we find no significant difference between 
the environments of the host galaxies and the environ- 
ments of comparison galaxy samples. In the following 
subsections, we examine how these results compare to 
the results of related studies in the literature. We also 
investigate potential selection effects, which could bias 
our supernova and galaxy samples, and finally we discuss 
the implications of our results in terms of the currently- 
unknown SN la progenitor population. 

4.1. Comparison to Previous Work 

Recent analyses have reached a variety of conclusions 
regarding the dependence of the SN la rate on envi- 
ronment. Using angul ar cross-correlation te chniques on 
data from the SNLS, ICarlberg et"al] HqoI) found that 
supernova hosts at 0.2 < z < 0.9 are more strongly clus- 
tered than a sample of field galaxies selected to have 
the same redshift and (i-band) brightness distributions 
as the hosts. However, inaccurate photometric red- 
shifts could, by overbroadening the measured galaxy red- 
shift distribution, dilute the measured angular-clustering 
strength for the galaxy sample relative to that of the 
supernovae, for which spectroscopic redshifts were ob- 
tained. Furthermore, the comparison sample is not se- 
lected to match the stellar mass, rest-frame color, or 
SFR distributions of the host galaxies. Thus, this result 
could also be attributed to the dependence of the SN la 
rate on galaxy properties, where the SN la rate increases 
with stellar mass and star-formation rate, and to the 

density and low— density environments have the same luminosity 
distribution. If type la events are more or less luminous in lower- 
density regions, then the inferred value of R would decrease or 
increase, accordingly. 
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Fig. 7. — We plot the distribution of overdensities for the blue (left) and red (right) galaxies in the SN-B and SDSS-B samples. All of 
the histograms are normalized to have equal area. The total number of objects per bin for the blue or red SDSS-B and SN-B samples is 
shown by the left and right axis scales, respectively. We find that the environment distribution for the blue type la SN hosts in the SN-B 
sample is skewed to lower overdensities relative to a sample of blue galaxies with similar luminosity, color, and redshift distributions. 



TABLE 2 

Summary of Mean Overdensities 



Sample 


<logio(l + <55) > 


'^<logio(l-l-i5)> 


Sample 


< logioCl + ^s) > 


0'<logio(l + <55)> 








SDSS-A 


0.070 


0.005 


SN-A 


-0.132 


0.059 


SDSS-A [blue] 


-0.073 


0.007 








SDSS-A [red] 


0.160 


0.006 


SN-B 


-0.023 


0.098 


SDSS-B 


0.079 


0.007 


SN-B [blue] 


-0.345 


0.053 


SDSS-B [blue] 


-0.083 


0.009 


SN-B [red] 


0.235 


0.153 


SDSS-B [red] 


0.196 


0.009 


SN-(J 


0.011 


0.094 


9DS9-C 


0.100 


0.009 


SN-C [blue] 


-0.338 


0.064 


SDSS-C [blue] 


-0.077 


0.012 


SN-C [red] 


0.248 


0.138 


SDSS-C [red] 


0.219 


0.012 


SN-D 


0.009 


0.101 


9DS9-D 


0.082 


0.007 


SN-D [blue] 


-0.171 


0.130 


SDSS-D [blue] 


-0.090 


0.009 


SN-D [red] 


0.174 


0.146 


SDSS-D [red] 


0.229 


0.010 


SN-E 


-0.017 


0.096 


9DS9-E 


0.087 


0.007 


SN-E [blue] 


-0.317 


0.058 


SDSS-E [blue] 


-0.074 


0.009 


SN-E [red] 


0.235 


0.153 


SDSS-E [red] 


0.216 


0.010 








SDSS-F 


0.077 


0.007 








SDSS-F [blue] 


-0.062 


0.009 








SDSS-F [red] 


0.170 


0.009 



Note. — The mean and the error on the mean of the overdensity distributions for the various SNe and 
galaxy samples. 



dependence of galaxy properties on environment, where 
more massive galaxies favor overdense environs.^ In fact, 
when weighting their fie ld galaxy sample by s tellar mass 
and star-formation rate. lCarlberg et al.l(|2008f) show that 
the clustering of the type la supernova hosts and field 
galaxies in the SNLS are in good agreem ent. A parallel 
analys is of data drawn from the SNLS by I Graham et al.] 
()2008l ). focusing on SNe la identified within galaxy clus- 
ters at intermediate redshift, found no significant dif- 
ference between the SN la rate in clusters from that in 
field ellipticals. However, their supernova sample is quite 
small (only three probable cluster la events) and the re- 
sults are dominated by statistical uncertainties. 

Studying su pernovae in the local (z < 0.2) Universe, 
I Sharon et al.l (2007) found the type la rate in nearby 
clusters to be roughly consistent with estimates of the 

^ There is also evidence that the most strongly star-forming 

S alaxies tend to reside in dense environments at 2 ^ 1 
Coo£er_e^^l. 2008b), which might play a role in biasing the 
ICarlberg et al. ( 2008,) sample, since it is strongly weighted towards 
SNe at 2 > 0.5. ' 



SN la rate in local elliptical galaxies; this work, however, 
was based on a sample of only six cluster supernovae 
and thus yielded very large uncertainties in the measured 
cluster type la rate. A more recent analysis of dat a in 
nearby {z < 0.04) clusters bv lMannucci et all ()2008D ar- 
rived at a very different result, finding that the SN la 
rate (per unit mass) is more than 3 times higher in clus- 
ter ellipticals relative to field ellipticals, using a sample 
of 11 cluster and 5 field type la events; an identical SN 
rate in b oth samples is exclude d with P ~ 0.02. In their 
analysis. iMannucci et al.l (|2008') compare the rate within 
galaxies of like mass and morphology, attempting to re- 
move the known correlations between galaxy properties 
and environment. For this reason, and the simila rity in 
redshift range probed, the Man nucci et al.l ()2008f ) study 
provides the most significant and relevant comparison to 
our results. 

In particular , the environment-dependent rates of 
llVlannucci et al.l Hool) are most closely connected to our 
results regarding the env ironments of red SN la hosts. In 
contrast to the results of llVIannucci et ahl (|2008f ) , we find 



The Environments of Local SNe la 



11 



TABLE 3 
Summary of Statistical Tests 



Samples 




-PwMW 


Pks 


SN-A/SDSS-A 




5.5 X 10-'' 


6.8 X 10-5 


SN-A/SDSS-A 


blue] 


0.011 


0.030 


SN-B/SDSS-B 


bluej 


0.005 


0.006 


SN-B/SDSS-B 


red] 


0.499 


0.846 


SN-C/SDSS-C 


blueJ 


0.006 


0.024 


SN-C/SDSS-C 


red] 


0.479 


0.632 


SN-D/SD9S-D 


blue] 


0.068 


0.057 


SN-D/SDSS-D 


red] 


0.336 


0.639 


SN-E/SDSS-E 


blue] 


0.011 


0.018 


SN-E/SDSS-E 


red] 


0.479 


0.850 


SN-B/SDSS-F 


blue] 


0.002 


0.005 


SN-B/SDSS-F 


red] 


0.435 


0.811 



Note. — We tabulate the P-values, PwMW 
and -Pks I from comparing the environment val- 
ues in the listed samples, using the one-sided 
Wilcoxon-Mann-Whitney (WMW) U test and the 
two-sided Kolmogorov— Smirnov (KS) test. As dis- 
cussed in smaller values indicate a lower proba- 
bility that the observed differences in the samples 
will occur by chance if they are selected from the 
same underlying parent distribution. Note that, 
as a one— sided test, PwMW has a maximum value 
of 0.5. 

no significant trend sucli that supernovae occur more of- 
ten in overdense regions; within our SDSS samples, the 
environment distribution for the subset of red hosts is 
indistinguishable (within the uncertainties) from that of 
the samples of red comparison galaxies (see Table [5] and 
Table [21). However, our results regarding the environ- 
ments of red host galaxies have large statistical uncer- 
tainties associated with them, making any strong state- 
ment about the environment dependence of the type la 
rate in red (or elliptical) hosts impossible. 

To investigate an y potential discrepancy between our 
results and those of iMannucci et al] ()2008l ). we attempt 
to test the likel i hood that we would have detected the 
IMannucci et al.l (|2008f ) result of a > 3x higher SN la 
rate in cluster ellipticals relative to field ellipticals, given 
our sample size of 25 red SN la hosts in SN-B. From 
the SDSS-B galaxy sample, we select the 4256 random 
red galaxies which match the rest-frame g — r color, ab- 
solute r-band magnitude, and redshift distributions of 
the red SN la hosts in sample SN-B. From this par- 
ent population, we then draw 5000 independent sam- 
ples of 25 random galaxies with the likelihood of drawing 
an object in a high-density (log]^g(l -I- S^) > 1) environ 
forced to be 3 x greater than that of a galaxy in a low- 
density (log]^Q(l + ^5) < 1) environ. For each realization, 
we compare the distribution of environments for the 25 
mock SN hosts to the full sample of red galaxies in the 
SDSS-B sample, using the Wilcoxon-Mann-Whitney U 
and Kolmogorov-Smirnov statistical tests. Finally, this 
mock sample-selection analysis is repeated with galaxies 
in high-density environments selected at 1 x , 2 x , and 4 x 
the rate of galaxies in low-density regions. 

As shown in Figure [3 in only ^ 50% of the mock re- 
alizations with a 3x higher likelihood of SNe occurring 
in dense environs, would we have rejected the hypothe- 
sis that the environments of SN hosts are on average no 
more overdense than that of the SDSS-B red galaxy sam- 
ple with P < 0.05. The likelihood that we would have 
detected a difference in the typical environment becomes 
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Fig. 8. — For 5000 mock realizations of a SN la rate 3x higher 
in high-density than in low-density environments, we plot the 
percentage of realizations that yield a one-sided Wilcoxon-Mann- 
Whitney (WMW) probability less than Px (solid black line) and 
a two-sided Kolmogorov-Smirnov (KS) probability less than Px 
(solid red line). The black dotted, dashed, and dashed-dotted 
lines correspond to mock realizations simulating a SN la rate in 
high-density regions 4x, 2x, and Ix that in low-density environs, 
respectively. Note that the black lines reach a fraction of unity 
at a probability (Px) of 0.5 since they correspond to the one-sided 
WMW test, while the red lines saturate at P^ = 1 since they corre- 
spond to the two-sided KS test. The black and red vertical dotted 
lines give the WMW and KS probabilities, respectively, computed 
from a comparison of the environments of red galaxies in the SN-B 
and SDSS-B samples (see Table [3)1. 

even smaller when we assume a SN la rate that is only 
2x higher in dense environs; in this case, the WMW test 
P-value is <0.05 in 25% of the mock realizations. 

When we apply the two-sided KS test, we find similar 
results; for the majority (> 60% and > 85%) of mock real- 
izations, we would not have detected a significant differ- 
ence in the environment distributions at a P > 0.05 sig- 
nificance level, assuming a SN la rate 3x and 2x higher 
in overdense regions, respectively. Thus, while our anal- 
ysis of type la supernovae in red SDSS ho s ts doe s not 
support the conclusions of IMannucci et all ()2008D . our 
results are not in direct conflict either, given the large 
uncertainties in both results. With that said, the high 
WMW and KS P-values computed from a comparison of 
the environments for the red SN-B and SDSS-B galaxies 
(see Table [3]) occur only rarely in our simulated galaxy 
samples^° and thus favor a SN la rate in cluster ellipticals 
more in line with that of fleld ellipticals, at the low end 
of the possible ran ge given b y Mann ucci et al.l (|2008f ). 

While the Mannucci et al.i (2008) work has a sample 
size roughly double that of iSharon et al.l (|2007[ ). the un- 
certainties in their measured SN la rates are still quite 
large, such that the field and cluster rates are only dif- 
ferent at a 2.1(7 level, thereby leaving much uncertainty 
in the dependence of the type la rate on environment 
within early-type galaxies. 

4.2. Potential Selection Effects 

10 We find PwMW = 0.499 and Pks = 0.846 in ~0.05% and < 
5% of the mock realizations with a 3x higher SN rate, respectively. 
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While the SDSS-II Supernova Survey includes > 100 
SNe la with spectroscopic follow-up in the redshift range 
0.05 < z < 0.15, less than half of the total sample 
is successfully matched to a host galaxy in the SDSS 
DR6 spectroscopic sample. Given the large number of 
unmatched supernovae, it is important to understand 
any potential biases that could arise from the particular 
method adopted to identify host galaxies. For example, 
our observation that SNe la in blue galaxies are biased 
towards low-density environments could naturally arise 
from employing an algorithm to select host galaxies that 
is for some reason biased against finding hosts in over- 
dense regions (e.g., groups or clusters). Given the strong 
correlations between environment and galaxy properties, 
it is also critical to understand any aspect of the host- 
identification methodology that could be biased towards 
identifying hosts of a particular galaxy color, luminosity, 
etc. 

To test the sensitivity of our results to the particular- 
ities of the host-identification algorithm, we define two 
additional SN samples, SN-D and SN-E (see Table [T] 
for sample definitions). The SN-D sample is selected to 
test for any bias in the SN-B sample towards SNe be- 
ing preferentially matched to red host galaxies. For the 
vast majority of SNe in the SN-B sample, only one pos- 
sible host galaxy is found within the cylindrical search 
window (see H2.3I for the details of the SN-B matching 
algorithm). However, for three SNe, multiple host galax- 
ies are identified within the SDSS spectroscopic data set. 
In Figure [TJ we mark the location of these possible, al- 
ternate ("secondary") host galaxies using magenta stars, 
with dotted lines connecting them to the location of the 
"primary" host, as defined in the SN-B sample. 

While the choice of the host galaxy in these ambigu- 
ous cases does not significantly impact the measured en- 
vironment of the SNe, the location of the host in color- 
magnitude space directly affects whether the host is clas- 
sified as "blue" or "red" in our color cuts and affects the 
composition of the comparison galaxy sample (SDSS-B). 
Confusion among hosts is likely to be more common in 
ovcrdcnsc regions, which could cause a small bias in the 
measured distribution of supernova environments when 
divided into blue and red subsamples. 

To test the robustness of our results to the ambiguity 
in the host identification, we apply a Murphy's Law ap- 
proach, where if anything can go wrong to bias our SN-B 
sample, then it will. In the ambiguous cases, we corre- 
spondingly match the SNe to the bluest host according 
to rest-frame g — r color. Given the potential host galax- 
ies, this creates the sample of SN hosts (SN-D) with the 
highest possible mean density among blue host galax- 
ies (see Table [J) . Even in this extreme (and somewhat 
unlikely) scenario, we still find that the distribution of 
environments for the blue SN hosts in sample SN-D is 
distinct from that of a corresponding sample of galaxies 
chosen to match in color, luminosity, and redshift (sam- 
ple SDSS-D), with the difference significant at a >94%, 
following a two-sided KS test. 

Another aspect of the host-identification methodology 

Two of the three SNe in SN-B with an ambiguous host iden- 
tification reside in overdense environments (log]^Q(l -|- 65) 1.5), 
and the third SN resides in an environment that is still more over- 
dense (logiQ(l + (55) ~ 0.36) than any of the SNe in the blue SN-B 
subsample. 



that might lead to a bias against identifying hosts for 
SNe in dense regions is the size of the cylindrical aperture 
used to search for possible host galaxies. If this aperture 
is too small, then the sample could be biased against in- 
cluding SNe in galaxies of greater physical size. Massive 
galaxies, which are inclined to resid e in overdense envi- 
ronments fe .g.. ICooper et al]|2008a| ). also tend to have 



larger sizes (IShen et al 



sions (iFaber fc Jacksoii 



2003 ) and larger velocity d isper 



19761 : iDjorgovski fc Davig [T987). 



which potentially could lead to SNe occurring at pro- 
jected and velocity separations outside of the windows 
used to define the hosts in SN-B. 

Recognizing these correlations between galaxy size, 
velocity dispersion, and environment, we alternatively 
identify a host galaxy sample (SN-E) using a projected, 
radial window of 40 kpc (physical) to identify poten- 
tial hosts on the plane of the sky in conjunction with a 
velocity window of Av = 5000 km s~^ along the line of 
sight, both ~60% larger than the windows of 25/i~^ kpc 
and Av — 3000 km s~^ used to define the hosts in sam- 
ple SN-B. Like when defining sample SN-B, when mul- 
tiple potential hosts are identified within this window, 
the galaxy closest in projected distance is taken as the 
host. Using these larger radial and line-of-sight search 
windows, only one additional supernovae is matched to a 
host (relative to the number of hosts identified in SN-B); 
however, the percentage of SNe unambiguously matched 
to a host declines slightly from 92% to 88%. Thus, we 
conclude that our SN-B sample is not significantly biased 
against SNe occurring in more massive galaxies. 

A final potential selection effect associated with the 
identification of SN hosts is the possibility that match- 
ing to the SDSS spectroscopic galaxy sample is biased 
against matching SNe to hosts in rich environments due 
to fiber collisions. Our SN-A and SN-C samples are 
specifically designed to test for such an effect. By match- 
ing directly to the imaging catalog, SN-C avoids any bias 
associated with the allocation of fibers in the spectro- 
scopic component of the SDSS. We find that the general 
results obtained by analyzing the SN-B sample are sup- 
ported by the analysis of the SN-C sample. Thus, we 
conclude that our results are likely robust to any bias 
associated with all ocation of SDSS fiber s in overdense re- 
gions. As shown bv lCooper et al.l (|2005( l and lGerke et al.l 
(|2005f) . dense regions on the sky do not always translate 
into dense regions in redshift space. 

Looking beyond selection effects that might bias the 
host-identification process, we also examine the possibil- 
ity that our supernova samples are biased against bulge- 
dominated systems. Given the higher surface brightness 
of bulges, relative to galactic disks, it can be more diffi- 
cult to obtain spectroscopic follow-up data for a super- 
nova near the bulge of a galaxy. For galaxies with larger 
bulges, this problem is obviously greater. Moreover, ob- 
servations of local type la SNe find correlations between 
the light-curve decline rate and the peak luminosity with 
morphology, such that late -type galaxies ho st brighter 
SNe la (e.g., iHamuv et aP [1995. i99& Gallagher et all 
l2005h . Thus, given the correlation between morphology 
and environment, where the bulge-dominated fraction 
increases with local density, our observed deficit of su- 
pernovae in high-density environments within the star- 
forming population and lack of an increase in the type 
la rate within cluster red-sequence members could both 
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be attributed to such a morphology bias. 

To test for this selection effect, we compare the mor- 
phology distribution of our host galaxies in the SN-B 
sample to that of the galaxies in the SDSS-B comparison 
sample. As a tracer of morphology, we utiliz e the Sersic 
indices as mea sured for SDSS galaxies by iBlanton et al.l 
(|2003bl l2005al) . While the Sersic index is a measure of 
morphology derived from the fit of only a single compo- 
nent to the galaxy's radial profile (e.g., versus bulge-disc 
decomposition), we find no significant difference in the 
morphologies of our supernova hosts relative to the mor- 
phologies of the comparison galaxies. 

In much the same way that it can be difficult to detect 
a supernova superimposed on a bulge relative to a disk, 
follow-up spectroscopy of candidate events (used to con- 
firm the supernova type) is also affected by the relative 
brightness of the host galaxy; spectroscopic observations 
of supernovae in brighter hosts are more difficult due to 
contamination of the supernova spectrum by emission 
from the host galaxy. Any incompleteness in our sam- 
ple that is dependent on apparent magnitude, however, 
is effectively controlled for in our samples by matching 
the comparison galaxy samples (e.g., SDSS-B) to the 
supernova host samples (e.g., SN-B) according to color, 
luminosity, and redshift. By matching both luminosity 
and redshift, the samples being compared also have the 
same brightness. For the very same reason, our analysis 
is insensitive to any incompleteness in the SDSS spec- 
troscopic catalog that may result from failing to obtain 
redshifts for fainter galaxies. 

Finally, we test the robustness of our results to the 
particularities of our adopted color division between blue 
and red populations. To do this, we shift the color cut 
given in Equation [1] by ±0.04 magnitudes in rest-frame 
g — r color, which results in changes in the blue and red 
components of the SN-B supernova sample of ±3 super- 
novae. With these changes in the subsamplc definitions, 
our results regarding the environments of blue host galax- 
ies remain "highly significant" (i.e., P < 0.01). Further- 
more, even when shifting the color cut by as much as 0.1 
magnitudes, the WMW U and KS tests indicate that the 
blue supernova host galaxies populate a significantly (i.e., 
P ^ 0.05) distinct distribution of environments from the 
blue comparison galaxies in the SDSS-B sample. 

4.3. Supernova la Progenitors 

As discussed in fJTl observational studies of nearby and 
distant supernovae have supported the definition of two 
components to the type la supernova rate, a "prompt" 
and a "delayed" component. However, there is no current 
observational evidence that indicates two distinct pro- 
genitor channels associated with the two components of 
the type la rate. With that said, recent theoretical anal- 
yses have found difficulty reconciling observations with 
models in which both the "prompt" and "delayed" la 
components result from single-degenerate events, where 
a carbon-oxygen WD accretes ma tter from a non- 
degenerate, compan i on star (e.g., lYungelson &: Liviol 
120001 iGreggiol 120051 : iPritchet et all I2008D . Allowing 
for the possibility of double-degenerate scenarios, in 
which two WDs are drawn together via angular mo- 
mentum losses r esulting from gravitational radiation 
(|Iben fc Tutukovlfr984: Webbink,1984 ) , some success has 
been found at predicting observed SN rates as well as 



the chemical enrichment of local galaxies (|Greggicil l2005t 
iMatteucci et al.ll2006l : fGreggio et al.ll2008f l. 

If the two components of the type la rate are somehow 
comprised of single- and double-degenerate events, then 
the "prompt" component is likely primarily driven by 
single-degenerate events, due to their shorter minimum 
timescale for occurrence. That is, single-degenerate 
events are favored for the "prompt" component, since 
the minimum timescale between formation of the pro- 
genitor star and occurrence of the supernova is on the 
order of the lifetime of a SMq star (i.e., ~ 30 Myr), 
while the minimum timescale for occurrence of a double- 
degen erate event is considerably longer 1 Gvr. lGreggiol 
[2005h . Thus, our results suggest that single-degenerate 
events are favored (or are more luminous) in low-density 
environs, given the assumptions detailed above. 

Now, while physical mechanisms such as galaxy merg- 
ers or harassment are more common in high-density re- 
gions such as groups and clusters, it is difficult to suggest 
a physical mechanism specific to low-density environs 
that would directly impact the evolution of stellar pop- 
ulations and influence the SN la rate or luminosity (e.g., 
by raising the binary fraction so as to produce more type 
la events). It is seemingly more likely that environment 
is strongly correlated with a galaxy property that traces 
differences in stellar populations, such as stellar or gas- 
phase metallicity or stellar age. In the following section 
f! j4.4p . we discuss this issue in more detail. 

4.4. The Role of Metallicity 

As discussed in detail in fJH we find that blue SN 
la host galaxies are only seen in comparatively low- 
density regions, suggesting that prompt supernovae la 
occur (or are found) preferentially in underdense envi- 
rons or that they are somehow suppressed in the high- 
density regime. However, there is little physical motiva- 
tion for connecting large-scale environment-specific pro- 
cesses (e.g., mergers, strangulation, etc.) with the gen- 
eration or suppression of type la supernovae. A poten- 
tially interesting, though currently poorly unconstrained, 
possibility could be that merger-induced star formation 
has somewhat different properties (e.g., a different initial 
mass function) than typical star formation (e.g., due to 
higher gas densities). The prevalence of mergers in envi- 
ronments such as galaxy groups could thus be associated 
with an environment dependence to the prompt la rate 
or luminosity. 

Another possibility could be the suppression of SNe 
within galaxies in clusters due to ram-pressure strip- 
ping, where the (metal-poor) outskirts of systems are 
preferentially confiscated by the IGM. The stripped 
gas and stars would then contribute to the intraclus- 
ter light a nd result in intracluster (i.e., intergalac tic) su- 
pernovae (|Gal-Yam et al.l[200llMaoz et al.ll200l . Such 
supernovae could be missed by narrow-field supernova 
searches that target i ndividual galaxies in clu sters rather 
than large fields fe.g.. lCappellaro et al.lll999f) . The imag- 
ing data from the SDSS-II Supernova Survey, however, 
covers a large and nearly continuous field, such that 
any intracluster SNe would be included in the supernova 
sample. Still, an intracluster SN would be less likely 
to be matched to a host galaxy, thereby mimicking a 
suppression of SNe within galaxies in overdense regions. 
Given the large uncertainties in the cluster SN la rate 
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and in the relat i ve co ntribution from intracluster events 
(jGal-Yam et alJl2003f ). the role of such supernovae re- 
mains largely unknown. With that said, the relatively 
low concentration of blue galaxies in local clusters helps 
to minimize the impact of intracluster events on our re- 
sults, as related to the preferential occurrence of SNe la 
in blue host galaxies in underdense environs. 

Alternatively, it could be the case that environment 
is correlated with a galaxy property or properties, for 
which we did not control in our analysis. A correla- 
tion between a given physical property and environment 
would cause the overdensities about our SN host galax- 
ies to be skewed relative to the comparison sample, if 
supernova host galaxies are a (relatively) biased tracer 
of that galaxy property. One likely galaxy characteristic 
to consider is metallicity. 

The metal abundance of a galaxy is often quantified 
in two ways: [1] by measuring the metal abundance 
in the interstellar medium (the gas-phase metallicity) 
and [2] by assessing the amount of metals locked up 
in the stellar population (the stellar metallicity). The 
gas-phase metallicity is a product of the recent (< 1 
Gyr) accretion an d star-formation history of the galaxy 
(|Finlator fc Dav^ [2008). In contrast, the stellar metal- 
licity, which is commonly inferred from fits to stellar 
absorption features in optical spectra, traces the met- 
als locked up in the old stellar population. Thus, stellar 
metallicity is a tracer of the integrated star-formation 
history and a measure of the gas-phase metallicity in 
the galaxy when the bulk of the stars formed. Here, we 
investigate the potential role of both stellar and gas- 
phase metallicity in biasing our supernova samples in 
star-forming galaxies towards underdense environments. 

To study the potential role of stellar me tallicity in our 
result s, we employ the measurements of iGallazzi et al.l 
((2005h . which are based on model fits to spectral ab- 
sorption features in the SDSS DR4 spectra. The mea- 
surements are sensitive to the signal-to-noise ratio of 
the spectrum (see Table 1 of IGallazzi eral|[2005t ). which 
limits our ability to directly constrain the stellar metal- 
licitics of the SN hosts. More than half of the SNe in 
our sample are at z > 0.1, which means that the host 
galaxies tend to be relatively faint in the r-band (less 
than 25% of the SNe have an r-band magnitude < 16.5) 
and the derived metallicity values are highly uncertain. 
For this reason, we are unable to make any meaning- 
ful statements about the stellar metallicity values of the 
hosts relative to the comparison galaxy samples. 

However, we are able to study the relationship be- 
tween stellar metallicity, Z, and environment for the 
general SDSS galaxy population, which allows us to un- 
derstand how our results regarding the environments of 
type la hosts could be understood in terms of a stel- 
lar metallicity bias. In the Appendix, we investigate in 
detail the relationship between metallicity and environ- 
me nt for three magnitud e-limited samples drawn from 
the IGallazzi et al.l (2005) catalog. As shown in Figure 
1111 we find a significant correlation such that galaxies 
with more metal-rich stellar populations typically reside 
in more overdense environs. Furthermore, this general 
correlation persists when focusing on just those galax- 
ies that reside on the blue cloud (following Equation [T]) ; 
although, the trend is considerably weaker. 

A corresponding relationship between gas-phase 



metallicity and envir onment was recently published by 
ICooper et al.l ()2008al ): studying star-forming galaxies in 
the SDSS DR4, they found a significant correlation be- 
tween average gas-phase metallicity and local galaxy 
density, such that more metal-rich galaxies favor re- 
gions of higher overdensity. Along the blue cloud, this 
metallicity-density relation is comparable in strength 
to the well-known c olor-density relation. Moreover, 
I Cooper et al.l ()2008al ) show that metallicity has a rela- 
tionship with environment separate from that observed 
with color and luminosity (or with stellar mass). Gas- 
phase metallicity is somewhat unique in this regard, as 
other galaxy properties (e.g., surface brightness, Sersic 
index, or stellar mass) that are strongly correlated with 
overdensity do not show a relationship with environment 
separate from that observed with color an d luminosity 
(jBlanton et al.ll2005aHCooper et al.ll2008af) . 

In order for the results in !j3]to be due to either of these 
observed correlations between environment and metallic- 
ity (either stellar or gas-phase), the blue SN host galax- 
ies must be biased towards lower stellar or gas-phase 
metallicities than the comparison galaxy samples. As 
discussed above, however, we are unable to directly test 
for any potential bias in our supernova samples, due to 
the lack of metallicity information (both stellar and gas- 
phase) for the majority of the hosts. Given the low 
signal-to-noise of the SDSS spectra, reliable metallicity 
measures are not feasible. Without direct constraints on 
the metallicities of the SDSS host galaxies, our results 
regarding the environments of SNe la in star-forming 
galaxies could still (at least partially) be understood in 
terms of a metallicity effect if the luminosity or the rate 
of type la events depends on metallicity, such that intrin- 
sically brighter supernovae arise from more metal-poor 
progenitors or such that the type la rate is elevated in 
metal-poor host galaxies. 

How ever, recent observational work by iHowell et al.l 
(|2008l ). studying SNe la in a relatively large sample 
(> 100) of star-forming and quiescent host galaxies at 
z < 0.75, suggests that variation in gas-phase metal- 
licity only accounts for a small portion (< 10%) of the 
measured dispe rsion in type la lu minosities (see Figures 
5, 7, and 11 of lHowell et al.ll2008i'). A variety of previ- 
ous studies (e.g .. Hamuv et al.ll2000l : llvanov et al.l 120001 : 
iGallaeher et all 120051) . using both gas-phase and stel- 
lar metallicity estimates, also found no significant evi- 
dence for a correlation between metallicity and the prop- 
erties of SNe la. With that said, whether low-metallicity 
stars might yield more luminous type la events is still 
a relatively poorly constrained question. For exam- 
ple, the gas - phase metallicity estimates employed by 
IHowell et al] (|2008! ) are derived from using measure- 
ments of the hosts' stellar masses to estimate the oxygen 
abundances according to the median mas s-metallicity re- 
lation for nearby star-forming systems (jTremonti et al.l 
l2004f ). This method is applied even to those systems 
thought to be quiescent and are thu s not included in 
the analysis of iTremonti et all (|2004( l. Future observa- 
tions, yielding direct oxygen abundance measurements, 
are needed to better determine the relationship between 
metallicity and type la luminosity. 

While the luminosities of type la SNe show no sig- 
nificant correlation with metallicity, our observations of 
supernova environments could also be explained by an 
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increase in the SN la rate in metal-poor hosts. From 
an observational standpoint, the dependence of the type 
la rate on host metallicity remains unconstrained, with 
current supernova samples including precision metallic- 
ity measurements for only a small portion of the host 
population. Still, recent observational studies show that 
a considerable number of type la s upernovae have been 
found in low-met allicity galaxies (jStrolger et al.l l2002l : 
iPrieto et al.l 120081 ). On the other hand, observations of 
type la supernovae in nearby, star-forming galaxies with 
significantly enriched interstellar media are evidence that 
the SN la rate is not zero in the mctal-rich regime (e.g., 
IGallagher et"an[2005h . 

In addition, theoretical models of supernovae suggest 
that the type la rate depends significantly on metal- 
licity, such that the ra te is lower in galaxies with 
lower metallicities (e.e,.. | Tornambe fc Matteucci 119861 : 
iHachisu et al.lll996tlKobavashi et alTl99ai200a but also 
see lUmeda et al.lll999b[) . This theoretical metallicity ef- 
fect, however, works counter to that which would explain 
our observations of supernova and galaxy environments, 
making our results regarding the environments of type 
la SNe in star-forming systems even more remarkable. 
With that said, the theoretical models are generally con- 
strained on a limited basis, being forced to match obser- 
vations of chemical evolution in the solar neighborhood. 

In spite of the theoretical predictions regarding the de- 
pendence of the SN la rate on metallicity, our results 
appear to be most easily understood in terms of a gas- 
phase metallicity effect, where prompt SNe la preferen- 
tially arise from metal-poor progenitors. In particular, 
the sharp cut-off in the observed distribution of envi- 
ronments for SNe la in blue hosts could be the result 
of a sharp feature in the metallicity distribution of the 
prompt SNe la progenitor population, such that prompt 
type la events rarely result from the evolution of rela- 
tively metal-rich stars. 

This picture is supported (though circumstantially) by 
several key points. First, the bias in the galaxy envi- 
ronment distribution for SNe la in star-forming systems 
towards underdense regions is likely attributable to the 
prompt (versus delayed) component of the type la pop- 
ulation, since delayed SNe are seen in both star-forming 
and quiescent galaxies and we find no evidence for any 
dependence of the type la rate on environment within 
red, non-star -forming galaxies. Fu rthermore, weak ev- 
idence exists (jMannucci et al.ll2008l ) to suggest that the 
SN la rate is actually higher among elliptical galaxies in 
clusters versus the field, which would work counter to 
the trend we observe in star-forming hosts, thereby sug- 
gesting that delayed and prompt SNe la have opposing 
relationships with galaxy environment, such that prompt 
events, which dominate in the star-forming galaxy pop- 
ulation, favor low-density environs and delayed events, 
which comprise all of the SNe observed in quiescent sys- 
tems, are preferentially found in high-density environs. 

There are two significant reasons that gas-phase metal- 
licity (and not another galaxy property such as stellar 
metallicity or age) is likely driving our observational re- 
sults. As discussed above, the prompt SN la component 
is correlated with star formation and thus thought to re- 
sult from the evolution of more massive stars (perhaps 
with masses of > 8M0). For this reason, the metal- 
licity of the progenitor population is more likely to be 



connected to the gas-phase metallicity and not stellar 
metallicity. 

In addition, in our analysis, we compare the environ- 
ments of our SN host galaxies to samples of galaxies 
with matched rest-frame color, luminosity, stellar mass, 
and star-formation rate distributions. Thus, for our 
results to be connected to a particular physical 
property, then that property must have a rela- 
tionship with environment separate from that ob- 
served with color and luminosity (or stellar mass 
or star formation rate) within the star forming 
population. Gas phase metallicity is the only 
galaxy characteristic know n to have such a r ela- 
tionship with environment ([Cooper et al.|[200 8a'). 
As discussed in more detail in the Appendix, stellar 
metallicity shows no relationship with environment sepa- 
rate from that observed with color and luminosity along 
the blue cloud. Moreover, even though type la luminosi- 
ties may depend o n stellar age in some systems (e.g., 
early-type svstems. IGallagher et al.ll2008l ). we also find 
that luminosity-weighted mean stellar age exhibits no 
significant correlation with environment at fixed color 
and luminosity among the star-forming population (see 
Appendix). 

While our results appear to be the manifestation of a 
metallicity bias in our supernova host samples relative to 
the comparison galaxy samples, such that the hosts are 
more metal-poor and the type la rate (or luminosity) is 
higher at lower metallicity, the exact physical explana- 
tion for why the SN la rate (or luminosity) would be ele- 
vated at low metallicities remains unaddressed. One pos- 
sible explanation for an increase in the SN la rate would 
be variation in the stellar initial mass function (IMF) 
with metallicity. However, measurements of the IMF in 
nearby star-forming regions show that it does not de- 
pend on nietallicity down to very lo w masses (< O.IMq, 
iBatd 120051 lYasui et al.ll2006L l2008f ). In addition, there 
is no observational evidence to suggest that the IMF 
varies with galaxy environment; instead, the stellar IMF 
is generally thought to be universal (at least locally, e.g., 
Elmegrecn 1999; Kroupa 2007; Sclman & Melnick 2008). 

Alternatively, a decrease in the SN la rate with metal- 
licity could be attributed to metallicity-dependent vari- 
ation in the binary fraction. Within the standard model 
for type la supernovae, an increase in the binary frac- 
tion should lead to an increase in the SN rate. How- 
ever, there is no evidence to suggest that the binary 
fraction shows any variation with metallicity (or environ- 
ment) across a broad range of me tallicities in the local 
Univ e rse (e.g., [C arnev et al."2005', but see also iMachidal 
120081 iMachida et al. 2009). Furthermore, type la su- 
pernovae might not result solely from the evolut ion of 
binary syst ems. As dis cussed in more detail by iMaod 
(|2008f ) and iToutI (|2005l ). our knowledge of type la su- 
pernova progenitors is quite poor, and a "single-star" 
SN la channel could help explain some observations of 
type la events that remain poorly understood within the 
standard binary-driven model. 

Finally, as highlighted throughout this paper, our re- 
sults are also consistent with an increase in the type la 

Morphology and surface brightness hkewise show no relation- 
ship with galaxy density se parate from that observed between color 
and luminosity l|Blanton et ai.ll2d05af ). 
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luminosity (rather than an increase in the rate) at low 
metallicities. While some early one-dimensional simu- 
lations of type la explosions suggested that there is a 
positive (or perhaps no) correlati on between metallic- 
ity and type la luminosity (e.g.. llJmeda et all Il999al : 
iDommguez et al.l 12001 ). more recent theoretical anal- 
ysis spanning a broader range of metallicities support 
a picture in which more metal-p oor progenitor s yield 
more luminous type la events (iTimrnes et al.l 120031 : 
iTownslev et"all l2009t iKasen et al.l I2009D. IiTparticu- 
lar, the analytical work of iTimmes et al.l (|2003[ ) con- 
cludes that the peak luminosity of type la events is 
anti-correlated with the metallicity of the pro genitor, 
such t hat the ^^Ni mass , a pro xy for luminosity (jArnettl 
Il982t iPinto fc Eastman! |2000D . produced in a type la 
event decreases linearly with increasing metallicity — 
i.e., luminosity depends roughly linearly on the gas- 
phase metallicity in the case of supernovae with short 
delay times, where the current gas-phase metallicity re- 
flects the stellar metallicity of the progenitor. So, while 
there is currently no observational evidence for a correla- 
tion between type la luminosity and gas-phase metallic- 
ity with in star-forming hos t galaxies (i.e., among prompt 
SNe la, iHowell et al.ll2008D . there is theoretical support 
to suggest that our results could be understood in terms 
of a luminosity effect in lieu of a bias in the type la rate. 

5. SUMMARY 

In this paper, we present an analysis of the large- 
scale environments of type la supernovae in the local 
(0.05 < z < 0.15) Universe, using data drawn from 
the SDSS-II Supernova Survey and from the SDSS DR6 
database. We estimate the local overdensity about each 
galaxy according to the projected fifth-nearest-neighbor 
surface density. For a range of supernova host galaxy 
subsamples, we define comparison galaxy subsamples, se- 
lected to match the color, luminosity, and redshift or stel- 
lar mass, star-formation rate, and redshift of the hosts. 
Using a variety of statistical tests, we then compare the 
distribution of environments for the supernova hosts and 
the comparison galaxies. Our principal results are as fol- 
lows. 

1. We find that type la supernovae in blue (i.e., star- 
forming) host galaxies tend to reside in low-density 
environments relative to galaxies of similar stel- 
lar mass and star-formation rate (or similar rest- 
frame color and luminosity). This analysis repre- 
sents the first constraint on the environment depen- 
dence of the SN la rate in star-forming systems. 

2. We find no difference between the environment 
distributions of red host galaxies and galaxies of 
like stellar mass and star-formation rate, within 
the observational uncertainties. This result is 
in some co ntrast to recent observational work by 
iMannuccTet al.l (j2008), which found a SN la rate 
more than 3x higher in local cluster ellipticals rel- 
ative to field ellipticals, though within the errors 
of each measurement. While our sample sizes are 
still somewhat small (;<30 la events in red hosts), 
tests on simulated galaxy catalogs suggest that the 
distribution of environments for red, SN la hosts 
presented herein is in poor agreement with a clus- 



ter SN la r ate as strongly e l evated relative to the 
field rate as iMannucci et all (|2008f ) find. 

3. We find a strong cut-off, at logjQ(l + (55) ~ 0.15, in 
the observed distribution of environments for type 
la supernovae in star-forming galaxies. Statisti- 
cal tests show that such a strong cut-off in over- 
densities is likely to result from a type la SN rate 
that is > 2 X higher within star-forming galaxies in 
low-density environments relative to like systems 
in more overdense regions. 

4. We conclude that variation in gas-phase metallicity 
is the most likely explanation for the observed dif- 
ference between the environments of type la SNe in 
star-forming galaxies and galaxies with like color, 
luminosity, stellar mass, and star-formation rate; 
type la events are more luminous or more numer- 
ous in metal-poor galaxies. Given existing observa- 
tional constraints, which indicate a lack of type la 
luminosity dependence on metallicity, our analysis 
suggests that prompt type la events preferentially 
result from the evolution of relatively metal-poor 
stars (i.e., have a higher rate in metal-poor sys- 
tems). 

5. Theoretical models of exploding white dwarfs likely 
need to be revisited, since they predict a strong 
fall off in the type la rate at low metallicity, while 
our observations unmistakably suggest the opposite 
effect for the prompt SN la channel. 

6. We predict an increase in prompt type la super- 
novae at higher redshift, as star-formation activity 
increases and metallicity decreases. 
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APPENDIX 

To study the relationship between stellar metaUicity , and environment for the general SDSS galaxy population, we 
select three subsamplcs from the Gal lazzi et al.l ()2005[ ) catalog down to limiting magnitudes of r = 16.5, 17, 17.5, which 
test for the dependence of our results on the precision of the metaUicity measurements while also probing increasingly 
more intrinsically faint samples. Figure [9] shows the distribution of the three magnitude-limited samples in g — r versus 
Mr color-magnitude space. Extending to fainter samples includes more metal-poor systems, especially at the faint, 
blue end of the blue cloud (see Figure fTO|) . 




-18 -19 -20 -21 -22 -23 



Fig. 9. — The rest— frame g—r color and r- band absolute magnitude distributions for magnitude-limited samples drawn from lGallazzi et al.l 
(gOOD. 

For these three magnitude-limited samples, we find a significant correlation between stellar metaUicity and envi- 
ronment such that galaxies with more metal-rich stellar populations typically reside in more overdense environs, with 
the trend relatively independent of the magnitude limit of the sample (see Figure fTT|) . This general correlation per- 
sists when focusing on just those galaxies that reside on the blue cloud (following Equation [T]) ; however, the trend is 
considerably weaker. 

As stated in H4.4[ for the main observational results of this paper to be connected to a particular physical property, 
that property must have a relationship with environment separate from that observed between color and luminosity 
(or stellar mass and star-formation rate). To test whether stellar metaUicity has such a relationship with environment, 
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Fig. 10. — The mean stellar metallicity, Z, as a fu nction of rest-frame g — r color and r-band absolute magnitude for galaxies brighter 
than r = 17.5 in the catalog of[Gallazzi ct al] 120051) . The mean metallicity is computed in a sliding box of width AM^ = 0.2 and height 
A[g — r) = 0.05, as shown in the upper left corner. 

we remove the dependence of environment on g — r rest-frame color and r-band luminosity for a magnitude-limited 
(r < 17) sample of galaxies on the blue cloud, using the methodology detailed bv lCooper et al.1 ((2008a). In short, we 
subtract the mean overdensity at the color and luminosity of each galaxy from the measured overdensity: 

As ^logioll + '^s)- < logio(l + '55)[g-r,A/,] >, (1) 

where the distribution of mean environment with color and absolute r-band magitude, < log]^o(l + ^5)[ff ~ '"i ^r] >, is 
median smoothed on A{g — r) = 0.15 and AA/r — 0.6 scales prior to subtraction. This difference gives the "residual" 
environ ment, A5, which quant ifies the overdensity about a galaxy relative to galaxies of like color and luminosity. As 
shown in lCooper et al.l (|2008af ). this new measurement of environment shows no correlation with color and luminosity 
(by design) or stellar mass and star-formation rate. 

We then study the relationship between residual environment and stellar metallicity, for our sample of star-forming 
galaxies with r < 17. As shown in Figure 1121 we find no significant residual trend between stellar metallicity and 
environment. That is, the mean residual environment is independent of stellar metallicity over the full range of 
metallicities probed. Thus, like morphology and surface brightness, stellar metallicity does not have a relationship 
with environment separate from that observed with color and luminosity. 

For the same sample of blue galaxies (with r < 17), we also examine the correlation between luminosity-weighted 
mean stellar age and environment. Again, we use the age measurements of Gallazzi et al. (200 5), which a re based on 
model fits to the spectra features in the SPSS spectra. The model spectra utilized bv lGallazzi et al.i (|2005f ) are derived 
from the population synthesis models of ,Bruzual fc CharlotI (|2003D and span a broad range of star-formation histories. 
The models are simultaneously fit to a minimum set of metal- and age-sensitivie spectral indices (D4000, H/3, H(5yi + 
H7A, [Mg2Fe], and [MgFe]'), yielding measurements of age and metallicity with typcial uncertainties on the order of 
at < 0.15 dex and az < 0.3. 

As shown in Figure I13i we find that there is a significant correlation between stellar age and absolute local envi- 
ronment along the blue cloud; galaxies with older stellar populations tend to reside in overdense environs relative to 
their counterparts with younger stellar composition. This effect is far from surprising given the correlation between 
stellar age and mass, wh ere more massive systems tend to have more enriched and older stellar populations (see Fig. 
8 of Gall azzi et aT][2005h . 

Now, we remove the mean dependence of environment on luminosity and color, again following Equation [U and 
study the correlation between residual environment, A5, and stellar age. We find that there is no sigificant relationship 
between age and residual environment (see Figure [M]) . That is, stellar age — like stellar metallicity — does not have 
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Fig. 11. — The mean st ellar metallicity, < Z >, as a function of local galaxy overdensity, logjQ(l + 65), in magnitude-limited samples 
drawn from lGallazzi et al.l lUoOS). 

a relationship with environment separate from that observed with color and luminosity. 
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Fig. 12. — The dependence of mean residual environment, A5, on stellar metallicity for the blue galaxies in a magnitude-limited (r < 17) 
sample drawn from Gallazzi ct al. (2005). The residual environment is computed following the methods of [Cooper et al. (2008^. We find 
no relationship between stellar metallicity and environment beyond that observed between color, luminosity, and environment. 




Fig. 13. — The dependence of mean absolute environment, (ogio(l + (Js), on stellar age, (t), for the blue galaxies in a magnitude-limited 
(r < 17) sample drawn from Gallazzi et al. {2003)- We find a significant relationship between mean stell age and environment such that 
older (i.e., more massive) galaxies favor overdense regions. 
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Fig. 14. — The dependence of mean residual environment, A5, on stellar age for the blue galaxies in a magnitude— limited (r < 17) 
sample drawn from Gallazzi ct al. (2005). The residual environment is computed following the methods of lCooper et al.l (I2008al) . We find 
no relationship between stellar age and environment beyond that observed between color, luminosity, and environment. 



